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1. introduction 

Nonlocality is one of the most profound aspects of a quantum mechanical system and 
it is a fundamental resource for quantum information processing. Nonlocality has 
been studied commonly in the operational perspective based on Bell's inequalities for 
bipartite two-dimensional systems with dichotomic measurements. The extensions to 
arbitrary dimensional systems have been proposed 01121 El- Recently, Kaszlikowski et 
al. considered joint probabilities of two distant measurements and suggested how to 
compare the strength of nonlocality between different dimensional systems for different 
numbers of measurement outcomes. A maximally entangled pure system violates 
Bell's inequality but as enough white noise is added the system loses its nonlocahty. 
Kaszlikowski et al. proposed the fraction of white noise as a measure of nonlocality, 
which may be used for its cross-dimensional comparison. Their numerical analysis for a 
maximally entangled state showed that the degree of violation increases monotonically 
with respect to the number of outcomes. 

More recently, Collins et al. 3j developed a family of Bell's inequalities for an 
arbitrary finite number of measurement outcomes. The family of inequalities are in 
good agreement with Kaszlikowski et al.^s results in terms of their measure of nonlocality. 
These studies would imply that the critical fraction of white noise is a useful measure 
in comparing the amounts of nonlocality for the different dimensional systems. The 
measure of nonlocality based on the noise fraction is, on the other hand, criticized by 
Acin et al. j3] as they found that partially entangled systems can give the larger violation 
(or stronger resistance to noise) of the Bell's inequality than the maximally entangled 
state. The other approach for a substantial violation of local realism was introduced by 
van Dam et. al. i5j and they found that CHSH inequality is the strongest nonlocality 
test for a bipartite system in terms of the statistical strength. 

In the test of Bell's inequality, a set of unitary transformations play an important 
role in the violation of Bell's inequality because local measurement settings for each 
party are characterized by local unitary transformations. In earlier studies [21 12 E], 
the transformations are restricted to quantum Fourier transformations (QFT). For a d- 
dimensional system the most general unitary transformation forms the group SU((i). It 
is thus questionable whether the QFT is sufficient to fully reveal quantum nonlocality. 
It has been known that the QFT is sufficient for a maximally entangled system of 
d = 2 and 3 . However, it is still an open question for other dimensional systems and, 
more importantly, for an arbitrarily entangled system. This question is investigated in 
this paper. 

Quantum nonlocality for continuous-variable (CV) systems has been studied in 
various contexts. Bell argued that the original Einstein- Podolsky- Rosen (EPR) state [7] 
would not violate Bell-like inequalities since it has a positive-definite Wigner function 
and thus its correlation function with respect to position and momentum observables 
can be simulated by local hidden variables. On the other hand, introducing dichotomic 
measurements such as even or odd parities of the photon number and presence or 
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absence of photons, Banaszek and Wodkiewicz [HI E] showed the nonlocahty of the 
EPR state. The measurements follow displacement operations, that is, translations in 
the phase spaces of the modes. However, the scheme by Banaszek and Wodkiewicz did 
not give a maximal violation for the inequality. This motivated Chen et al. jTU] to 
investigate another type of observable, with the unitary transformations other than the 
displacement in phase space, which results in the maximal violation for EPR state. The 
observable, so-called "pseudo-spin" operator, is defined as tensor summations of Pauli 
spin operators, which is element of an SU(2) group. 

The generalization of a dichotomic measurement to an arbitrary finite number 
of outcomes for a nonlocality test of a CV state was proposed by Brukner et al. 
[T^ as establishing a correspondence between a CV and a discrete system of an 
arbitrary finite dimension. However, in their work, it is not clear whether the 
correspondence can be given as a physically plausible map, i.e. a completely positive 
(CP) map. Moreover, in their analysis for nonlocality, Brukner et al. did not employ 
general transformations in SU((i) but the simple QFT transformations in varying the 
configuration of measurements. Thus the question of QFT being sufficient to reveal 
nonlocality arises in CV systems as well as in the finite-dimensional systems. 

In this paper, we formulate the Bell's inequality in terms of a linear summation of 
correlation functions which utilize the most general projective (i-outcome measurements. 
For the correlation function, we introduce a general form of correlation observable 
between two ci-level systems and find the eigenvalues from the generic conditions that 
any correlation function should satisfy. For the observable, we exploit all the possible 
unitary transformations in the SU((i) group on the configurations of local c?-outcome 
measurements. The subgroup algebra of SU(c/) allows us to prove that {(f — d) number 
of real parameters are sufficient to describe the local unitary operation. By inspecting 
symmetries and performing numerical analysis, we show, while the QFT suffices for 
a maximally entangled system, a partially entangled system requires more general 
transformations in order to fully investigate its nonlocality. 

The Bell's inequalities are applied to a CV state whose infinite-dimensional Hilbert 
space is decomposed onto the tensor sum of d-dimensional subspaces. The decomposition 
maps any CV state onto an arbitrary finite-dimensional state. We prove that the 
mapping is linear, trace preserving, complete positive. After applying the mapping, we 
investigate the violation of Bell's inequality for the two-mode squeezed vacuum state, 
as an example of a CV state, with the different outcome measurements. To search for 
the optimal violation in SU((i) transformation, several numerical methods are assessed. 

2. Bell's inequalities with (i-outcome measurements 

In this section, we investigate Bell's inequalities by considering the SU((i) group of 
transformations for the measurement with d outcomes. The series of inequalities may 
be derived by introducing "classically correlated observables" which can be constructed 
by local measurements and classical communications. 
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2.1. Special unitary transformation for the d-outcome measurement 

A measurement of a system is represented by a Hermitian operator which is called 
an observable. Any Hermitian operator on a (i- dimensional Hilbert space Tid can be 
expanded by the identity operator and the group generators of SU((i) algebra. Such 
a typical description in terms of group generators was introduced by Hioe and Eberly 
[Ti] . In order to obtain the generators of the SU((i) group, one may introduce transition- 
projection operators 

Pjk=\j){kl (1) 

where is an orthonormal basis set on Tid- Now, the (rf^ — 1) Hermitian operators 
are constructed as 

Ujk = Pjk + Pkj (2) 

Vjk = i{,P'jk - Pkj) (3) 




=->/|(7^(gA.-/P«,«) (4) 

where 1 < I < d — 1 and 1 < j < k < d. It is easy to check that when d = 2 these 
generators are Pauli spin operators. 

The set of G = {^12, Mis, . . . , vu, Vis, • • • , Wi, . . . , Wd-i} is composed of generators 
for SU((i) group, fulfilling the relations of tracelessness Tr(sj) = and orthogonality 
Tr(sjSj) = 26ij for Si, Sj G G. The elements Si E G hold the algebraic relation, 

[sj,Sk\ = 2iJ2fjkiSi (5) 

where fjki is the antisymmetric structure constant of the SU((i) algebra. 

The set G can be divided into the three mutually exclusive subsets of operators: 
U = {ujk}, V = {vjk}, and W = {wi}, which contain d{d — l)/2, d{d— l)/2, and (d — 1) 
elements respectively. The operators in each subset satisfy the algebraic relations 

[uij, Uki] = - iiSji (1 - Sik) Vik + Sik (1 - Sji) Vji 

+ 6jk (1 - 6u) Vii - 6u (1 - 6kj) Vkj}, 
[vij,Vki] = - i{Sji{l - Sik)vik + Sik{l - 6ji)vji (6) 

- 6jk (1 - Vii + 5ii (1 - 5kj) Vkj}, 
[wi,Wj\ =0. 

The commutators between the elements from U 01 V result in the operators in the 
subset V while any elements in W commute to each other. The commutation relations 
among the operators from the different subsets can be found as 

[Uij,Wi] = -{Xii-Xji)Vij 

[vij,wi] = - i{xii - Xji)uij (7) 



\u 



Vki] = iiSjkUii + SikUji - SjiUik - SiiUkj) 
+ 2z6ikS,i{\j){j\-\t){i\) 
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where the coefficient xu is given by 

and 6 is the Kronecker delta function. As a summary, one may symbohcally express the 
commutation relation among the subsets U, V, and W as 

[U, U] oc V, [V, V] oc V, [W, W](xO 

[U, W] oc V, [V, W] oc U, [U, V]ocU + W. (9) 

Using the SU((i) group generators, any Hermitian operator on the d-dimensional 
Hilbert space is represented by 

f^(a) = ^i + (10) 

where oq = Tr(](a) and aj = TrsjCl{d) are real numbers due to the hermiticity condition 
of the observable Cl{a). The coefficient aj comprises a {(P — l)-dimensional vector 
a = (ai, . . . ,arf2_i) which we call a generalized Block vector, while oq is constant over 
any SU((i) transformations. 

In the Heisenberg picture, the unitary transformation U of the Hermitian operator, 

n{a) n{a) = UVL{a)U\ (11) 

can also be described as a transformation of the generalized Bloch vector a. 
Decomposing Vt{a!) in the form (jlOj) with coefficients a'^ and using the invariance of 
the trace under cyclic permutation, the components a'^ of the transformed generalized 
Bloch vector are found to be 



a 



'. = Tr(%fi(a')) = TT{U^sfih{a)). (12) 



Since tj^^sfi is also Hermitian and traceless, it can be expanded in terms of the SU((i) 
generators as 

U^sfi = Y.Tjkh (13) 

k 

where Tjk = ^Tt(U^ sjUsk) is an element of a (rf^ — 1) x (d^ — 1) real matrix. The 
matrix T represents the direct relation between the transformed and untransformed 
generalized Bloch vectors a'^ = J^jTjkaj. As the norm of the generalized Bloch vector 
remains constant under the transformation, the real matrix T is orthogonal. 

An operator U in SU((i) can be represented in terms of the group generators 
s = (si,S2, ■ ■ ■,Srf2_i) as 

U{p) = exp{—ip- s) (14) 

where pis a (d^ — l)-dimensional parameter vector. The parameterization in Eq. fll4j) is 
said to be canonical. Experimentally, for the optical device, it is possible to realize the 
discrete unitary operation in SU((i) using biased multiport beam splitters ^Hj- In order 
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to derive the exphcit matrix elements of T in Eq. (|13|) corresponding to the unitary 
operator U{p), one may consider a set of differential equations for the generators: 




(15) 

where Sj{t) = JJ'^{tp)Sjtj{tp). After solving the differential equation and setting t = 1, 
the matrix T is derived in terms of the parameter vector p and the antisymmetric 
structure constant fkji as 

T(p) = exp(-2F(p)), where Fji{p) = J^Pkfkji- (16) 

k 

The antisymmetric characteristics of the structure constant fkji is related with the 
orthogonality of T as T^T = TT^ = Id'^^i where ]lrf2_i is an identity matrix on {(P — 1)- 
dimensional vector space. 

It is notable that a commutation relation appears in Eq. ()15|1 . From the fact that 
the group generators {sj} are divided into three subsets U, V and W and the generators 
of each subset satisfy the algebraic relations (jB)) and (jZj), one can find some symmetries 
for the rotation of generators. Especially, since [wi, Wj] = 0, it is possible to find that the 
rotation of generators {wi} G W along the direction of any Wj results in the generator 
itself as 

exp ( -ipyj^ Wj ) wi exp {ip^^ Wj) = wi (1 7) 

where pwj is the Wj component of the parameter vector appeared in Eq. ()14|1 . 

Eq. fll7|) implies that the dimensionality of the nontrivial parameter vector p for 
the unitary transformation of the Hermitian operator (fTT|) can be reduced to {d"^ — d). 
Without loss of generality, the Hermitian operator (fTIH) can be written with a given 
orthogonal basis {|j)} and the unitary operator of the basis transformation as 

n{d) = ij{pa)jZ^,\3){j\ U\Pa) (18) 

where Vtj is the non-degenerate eigenvalue. The generators in the subset W are sufficient 
to reconstruct all the diagonal bases that is, 

\3){3\ = \i-T.9lw3~i+k (19) 

fc=0 



where gl, = {1 — j^ko)^ 2{j+k)(j+k-i) • With help of Eqs. (fT7|) and (fT^ . one can see that 
the dimensionality of the nontrivial parameter vector Pa in Eq. (jl8|) is [d^ — d). This 
implies that any unitary transformation in SU(d) for the observable in the d-dimensional 
Hilbert space is sufficient with {d"^ — d) number of real parameters. 
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2.2. Classically correlated observable with d-outcome measurement 



Two observers, say, Alice and Bob perform local measurements on their own d-level 
systems and they communicate their outcomes via a classical channel. A classically 
correlated observable is thus constructed as assigning a weight /ijj for the pair of 
outcomes, i and j: 



E 



(20) 



^^iA^)a{A ® \j)b{i\- 

The correlation coefficient matrix n is a d x d real matrix. It is notable that Eq. (j20|) is 
the most general form of a correlation measure between any two d-\evel systems. The 
correlation observable E involves d"^ correlation coefficients. We show that without loss 
of generality, the number of independent parameters reduces to d and we determine 
their values. 

We require that a correlation observable should satisfy the following conditions. 
C.l A correlation function should be indifferent to local polarization, which means that 
TrEpA 0tB = TtEIa ® = 0, (21) 
Tr^^^p are the reduced density operators. This raises the following 



where pa,b 
condition: 



fiij = 0, \/i and ^ = 0, Vj. 



(22) 



For the case of two outcomes, there is the well accepted correlation matrix fi = 
{{yUi,i,/^i,2}, {/^2,i5y"2,2}} = {{1> -l}) !}}• Here, the translational symmetry, 
fJ'1,1 = Ai2,2 and /ii2 = P21, and equal spacing, /ii^i — ^2,1 = 2 leads the correlation 
observable to optimize the measure of correlation. We generalize these in the following 
two conditions. 

C.2 The correlation coefficients are unbiased over their outcomes (translational 
symmetry within modulo d): 

l^i+k,j+k = V/c. (23) 

C.3 The coefficients are equally separated and normalized (maximal discrimination): 

2 



for z > j 



d-l 

The condition C.2 leads the correlation matrix /i to be in the form of 

/ /il /i2 /U3 



(24) 



A^i A^2 



/id \ 



\ fJ'3 At4 

and further the condition C.l implies that 



(25) 



(26) 
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The condition C.3 determines all the /i;'s such that /ii = 1 for a maximally correlated 
state, fid = —1 for a maximally anti-correlated state and the other /i;'s are assigned to 
have equally spaced values between 1 and -1. Thus, the three conditions C.l, C.2, and 
C.3 uniquely determine the correlation matrix /i, 

(i - j)mod d 

^^Lj = 1 - 2 . (27) 

Using Eq. (fT^ . the correlation observable E in Eq. (j^UI) can be written in terms of 
the SU((i) generators as 

d-l 

E = P'k^iWk ® wi (28) 

k,l=l 

where fik,i is the transformed correlation matrix from fij^k- That is, 
fc+l z+i 

h,i = EE 9l-i+i9i-j+if^i,j (29) 
i=i j=i 

where gl is given in Eq. p9|) . Note that the correlation observable E in Eq. ()28|) does not 
contain any local identity operator 1^, due to the condition C.l. Further, the observable 
transformed by local unitary operations is written as 

E{p, q) = U{p) ® U{q)EU\p) ® U\q) (30) 

= E hm{P^q)Sl ® Sm 
Z,m=l 

where jji{'P-,q) = {p)fiT{q). The unitary operators (j{p) and tf{q) determine the 
measurement configuration for each side. Without any constraint for the li-outcome 
measurement, the unitary operators are subjected to the SU((i) group. 

2.3. Bell's inequalities for bipartite d- dimensional system 

In order to investigate nonlocality of a bipartite system, we introduce a Bell function 
which can be constructed by a linear combination of correlation functions of two parties. 
The Bell function can be written without loss of generality as 

B = Yc,E{p,,qi) (31) 

i 

where the correlation function E{p, q) = TtE{p, q)p and c = {cj} is an arbitrary vector 
which satisfies a normalized condition I]j Cj = 2 to make the Bell function B a polytope 
[in]. Note that the correlation function E{p, q) E [—1, 1] for all p and q. The classically 
correlated observable E{p, q) can be written as 

E{p, q)=Y l^^Aiv) ® PM, (32) 

where the projector Pi{p) = U{'p)\i){i\U'^{'P) is for the i-th outcome with the 
measurement configuration p and the correlation matrix fi is given in Eq. (fTTj) . The 
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joint probability that Alice and Bob obtain the outcomes i and j with the measurement 
configurations p and q is given by 

n,{p, q) = Tr (A(p-) ® P,{q)p) ■ (33) 

This implies that, from the joint probabilities for a given measurement, one can obtain 
the correlation functions for different measurement configurations and thus the Bell 
function (jHT|) . 

A Bell function has its boundary which is allowed by a local realistic model. It 
is worthwhile mentioning that quantum-mechanically correlated states do not violate 
the boundaries of all the possible Bell functions (jS!))- Only the Bell functions whose 
boundaries are violated by quantum-mechanically correlated states are of interest in 
the test of nonlocality [T7j. In this paper, we do not try to find all the classical 
boundaries. Instead, we consider the Bell function whose classical upper bound is 2 
with the particular vector c = (1,1,1,— 1), 

B = E{Ai, Bi) + E{A2, B2) + E{B2, ii) - ^(^2, ^i). (34) 

After a little algebra, one realizes that the Bell function (jMj) is exactly the same as 
the Bell function of Collins-Gisin-Linden-Massar-Popescu (CGLMP) 0, whose classical 
bounds are found as 2 with help of joint probabilities. 

Note that the third term of the correlation function in Eq. (jH^ has the parameters 
for the measurement configurations exchanged. In general, the correlation function 
E{p, q) depends on exchanging the parameter vectors, 

E{p,q)^Eiq,p). (35) 

The correlation function is invariant for the parameter exchange only for dichotomic 
measurements in which case the Bell function ()34|) is led to the CHSH Bell's inequality 

m- 

In order to find the quantum mechanical maximum for the Bell function Eq. 
CGLMP used the QFT unitary transformation, 

UQFT{A) = ^Y.e''^'^'^'^'^\j){kl (36) 
vet J I, 

which has only a single parameter (p^ to be adjusted for the measurement configuration. 
They found that for the ci-dimensional maximally entangled state their Bell function, 
which is the same as Eq. (jS^), has its maximum 

'^^^ ( 21 \ 1 
^ ^"^S V 2#sin^[7r(/ + l/A)/d] ^^^^ 

when {(pAi, <pA2y (pBi, (PB2) = (0)1/25 1/4,-1/4). This is always larger than the local 
realistic upper bound 2 and increases as the number, d, of measurement outcomes 
increases. However, as a special subset of the unitary group U{d), it is unclear whether 
the QFT measurement is optimal for the test of nonlocality when d > 3 even though it 
has been known that this is the case for maximally entangled states of d = 2 and d = 3 
0. 
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The raised question becomes rather dramatic if a state is partially entangled. For 
example, when d = 2, the violation of the Bell's inequality in Eq. (jHT|) is plotted in 
Fig. ^ Note that the Bell's inequality becomes the CHSH Bell's inequality when d = 2. 
The state is assumed to be in a pure state of I?/') = cosv9|00) + siny^lll). The Bell 
functions are optimized for the different measurement configurations : the dashed line 
is obtained by the QFT and the solid line by the SU(2) transformations. The figure 
shows that the QFT is not an optimal transformation in revealing the nonlocality of the 
partially entangled state. It is required to consider the general SU(c/) transformations 
for the optimal nonlocality test when a state is in a partially entangled state. 




Figure 1. Bell function B with respect to the entanglement e — ^/2sm(p for a two- 
dimensional bipartite state. Solid line represents for the case of SU(2) measurement 
and dashed line for the measurement with QFT. 



3. Highly degenerate measurement for a CV state 

In this section, we consider a CV state as a system to test its nonlocality. For the 
purpose of the nonlocality test, one needs to introduce a proper measurement which can 
show the violation of any Bell's inequality. Generally, for the case of a CV state, the 
spectrum of its measurement is continuous and the number of non-degenerate eigenvalue 
is infinite. Therefore, a difficulty arises, even in principle, in measuring the infinite 
number of outcomes and the test of its nonlocality. Several possible methods which 
can overcome such a difficulty have been suggested UHl 112] • These methods have 
adopted measurements which give a finite number of outcomes from a CV state. This 
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measurement naturally assumes an infinite degeneracy in the measurement. Recently, it 
is also found that homodyne measurement after single photon substraction from a CV 
state can play an essential role for a loophole-free nonlocality test [E]- 

In this section, we formulate the explicit form of an observable A which can give 
a finite d-number of outcomes from the measurement on a CV system. The observable 
corresponds to a mapping from a CV state to an arbitrary d-dimensional system and 
the mapping is the same mapping which was suggested by Brukner et a/. P^. We show, 
here, that the mapping is a linear, trace preserving, and complete positive (CP) map 
which implies that the density matrix of a CV state can be legitimately transformed 
into a finite dimensional state. 



3.1. d-outcome measurement for a CV system 



An observable A{a) which gives d outcomes from the measurement on a CV state can 
be found as a direct sum of the infinite number of the d-dimensional observables; 



i(a) 










fi(a) 









{31 



J 



where Cl{a) is the observable which is in the d-dimensional Hilbert space and has the 
explicit form as was given in Eq. (|1(J|). The measurement with the observable A{a) 
produces infinite degeneracy in each outcome since it counts every d modulo basis state 
as the same outcome. Alternatively, the observable A{a) can be written as 



oo d—1 



A{d) 



J2 ^jkia)\dm + j) {dm + k\ 

m=0 jk=0 



(39) 



where fljk{a) is the matrix element of the ci-dimensional observable which is 
parameterized with the d"^ — 1 dimensional generalized Bloch vector a. 

With the observable A(a), one can establish the mapping between a CV state and 
an arbitrary finite-dimensional state for the CV state. It exploits the fact |T2j that, 
from the physical perspective, any two systems can be considered as equivalent, if the 
probabilities for outcomes of all possible future experiments performed on one and on 
the other are the same. Mathematically, the requirement can be expressed as 

Tr (i(a)p) =Tr(^](a)prf) (40) 

where p is the density matrix of any CV state while pd is that of a d-dimensional state. 

Moreover, it is important to clarify that the mapping is a physically possible 
quantum process. One can show that the observable A{a) on a CV state results in 
a trace preserving, linear, CP map e as 

P^Pd = e{p) (41) 

where p G BiTi) and pd G B{1C). We denote that BiTi) is the set of operators defined in 
Ti. Note also that Ti and K, are the infinite and d- dimensional Hilbert spaces respectively. 
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In order to prove it, it is possible to make use of the correspondence between the complete 
positive maps and positive-semidefinite operators [201 • The density matrix in i3(/C) can 
be expressed by the transformation e as follows 

= Tr^^ (Ik; ® p^4) (42) 

where is a positive-semidefinite operator defined in B{}C ® Ti.). The positive- 
semidefinite operator has the explicit form as 

oo d— 1 

Re = Y.J2\'^)i^\^\dr^ + k){dn + l\ (43) 

n=0 kl=0 

which satisfies the trace preserving properties of the CP map by Ttk^{R^) = Ifi. The 
correspondence between the CP map and the observable A{a) is confirmed from the 
dual map of the map e |21j on the observable Q{a) and it is 

e''{n{a)) = Trie {^(a) ® iTi^J^) 

= A{a) (44) 

where denotes partial transposition on the Hilbert space Ti only. We conclude that 
the measurement with the observable A{a) on a CV state is equivalent to consider 
the state as a ci-dimensional state, which is mapped from the CV state, with the 
measurement of Cl{a). It can also be said that the mapping is a linear, trace preserving 
CP map from Eq. 



3.2. Mapping of multi-mode state 

The mapping B(T-Ca ® Hb) — ^ B{}Ci ® /C2) for a two-mode CV density matrix onto a 
bipartite (i-dimensional state is possible as 

P12 = Ttab (l/c ® l;c ® pIbRaiRb2) . (45) 

The mapping for an arbitrary number of modes can also be found as an extension of 
Eq. (gg). 

As an example, we consider a two-mode squeezed state which can be generated 
by a non-degenerate optical parametric amplifier 



, ^(tanhr)'^, , 

where \n) is a Fock state and r is the squeezing parameter. It is well-known that when 
squeezing parameter r goes to infinity, the two-mode squeezed state approaches to the 
EPR state jTj. From Eq. (jl^j) . one can map the two-mode squeezed state onto the 
d-dimensional pure state: 

l^Pa) = I ^^^^^ = ^(tanhr)"|n,n)A,B. (47) 
yi - tanh^'^ r n=o 

The mapped state is a partially entangled pure state whose entanglement is characterized 
by the squeezing parameter r. The state becomes separable only when r = and it 
becomes maximally entangled for the limit of r ^ 00. 
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4. Numerical analysis based on SU((i) group 

We investigate tlie optimal violation of the Bell's inequality based on the Bell function 
in Eq. (jMj) for the two-mode squeezed state. In order to search for optimization 
values of the inequalities, we employ several numerical methods such as steepest descent, 
conjugate gradient, and dynamic relaxation. Each method has its own advantages and 
disadvantages depending on the situations for optimization. The conjugate gradient 
leads to rapid convergence for a nearly hyperbolic function (where a bounded function 
looks like near its minimum). The steepest descent method enables one to find 
persistently lower values, even though it has disadvantages of slow convergence for the 
nearly hyperbolic function that is squeezed in parameter space. The dynamic relaxation 
method is in between the two methods. We consider the dynamic relaxation method in 
detail as the algorithms and implementations of the other methods can easily be found 
in literatures [231 ■ 

The dynamic relaxation method simulates a physical system under a potential and 
a friction, which resembles the Car-Parrinelo method for ab initio molecular dynamics 
|24j . Consider a bounded function B{{pi}) in terms of the parameter vector pi. For an 
optimization the method simulates a dynamic equation for a fictitious classical particle, 
by regarding Pi as its trajectory vector and B{{pi}) as a potential. The dynamic equation 
of motion can be written as 



where m is a mass of the fictitious particle and 7 is a friction ratio. Note that the 
equation is a kind of the Langevin equation. The particle will relax to the minimum 
of the potential. The solution to Eq. ()48|) approaches to the minimum of the function 
B[{pi}) exponentially due to the friction 7. A minimum is claimed to be achieved when 
\^p^B{{pj})\ X \pi\ < 10^®. For the numerical implementation a Runge-Kutta method 
is used to solve the dynamic equation with the following ranges of the parameters: 
m = 0.1, 7 G (0.5, 1.5) and St G (0.01, 0.1). The maximum value of B{{pi}) is obtained 
by replacing the "potential" B{{pi}) with —B{{pi}). For the optimizations of Bell 
functions in Eq. (jH^ . the parameter vector is given by a = {Ai, A2, Bi, B2) where Ai 
and Bi are parameter vectors for unitary transformations U (Ai) and tJ{Bi), respectively, 
in the group SU((i). 

We optimize the value of the Bell function B^ in Eq. ()34|) with SU((i) 
transformations for the two-mode squeezed state. All the results are checked and they 
are reproduced by conjugate gradient, steepest decent, and dynamic relaxation methods. 
Fig. 121 presents the optimized value of the Bell functions B^ with respect to the strength 
of squeezing tanh(r) for different numbers of measurement outcome. The Bell function 
is upper bounded hj 2, B^ < 2, under the local realistic theory. 

In Fig. 121 we note that a two-mode squeezed state always violates the inequality 
-Bd < 2 for all r > regardless of the number of measurement outcomes. Brukner et 
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Figure 2. Violation of the Bell's inequality based on the Bell functioned for two-mode 
squeezed state. Finite number of measurement outcome are considered. 

Table 1. The largest optimum values, Bd{rrn), of the Bell fimction for two- mode 
squeezed state, is the squeezing parameter which maximized the value of Bd for 
the given number, d, of measurement outcome. Bd(oo) is the value of B4 for an infinite 
squeezing. 
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al. calculated the values of the Bell function Bd for d = 3 based on the QFT for a two- 
mode squeezed state. They do not always achieve > 2 for the squeezing parameter 
r > even though a two-mode squeezed state is inseparable. Thus in order to properly 
achieve the optimum value of the Bell function, we have to consider all the possible 
transformations in SU((i). 

On the other hand, in the limit of r — > 00, a two-mode squeezed state becomes 
a regularized EPR state which is mapped onto a maximally entangled state in finite 
dimensional Hilbert space. However, for the maximally entangled state, the QFT suffices 
to obtain the optimum vahic of B^ as we have already discussed. 

For a given number of measurement outcomes, the amount of violation increases 
and decreases with respect to the squeezing parameter r. Let rd denote the value of 
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squeezing parameter that gives the largest violation for a given measurement with d 
number of outcomes. As shown in Tabled except for the dichotomic measurement, the 
infinitely squeezed state violates the inequality less than some partially entangled states. 
As increasing the number of outcomes, the largest optimum values of monotonically 
increase. It is also found in Fig. El that for the high squeezing regime the higher number 
of outcomes gives stronger violation while the result is reversed for the small squeezing 
regime. 

5. Final remarks 

We studied the most general d-outcome measurement for the Bell's inequalities of a 
bipartite system. In order to construct the inequahties, we introduced a classically 
correlated observable which is constructed in terms of local measurements and classical 
communications. For the configuration of the local measurements, we considered general 
transformations in SU((i). It was found that the number of parameters for the nontrivial 
operation is reduced to (d^ — d). After inspection of symmetries, we derived the Bell 
function that is composed of the correlation functions. This Bell function is equivalent to 
that found by Collins-Gisin-Linden-Massar-Popescu 0. The present numerical analysis 
shows that, when the system is in a maximally entangled state, the QFT is an optimal 
transformation for each local measurement. However, we show that this does not hold 
when the system is in a partially entangled state. 

In order to utilize the CV state for the nonlocality test, we investigated the mapping 
between a CV state and arbitrary dimensional system which was devised by Brukner 
et a/. |12. . We found the mapping is linear, trace preserving and complete positive map 
and it corresponds to a highly degenerate d-outcome measurement on a CV state. By 
applying the highly degenerate measurements, we investigated the optimal violation 
of the Bell's inequality for the two-mode squeezed state. Regardless of the degree of 
squeezing and the number of outcomes, the two mode squeezed state always violates 
the Bell's inequalities. This opens a possibility to extend Gisin's theorem jT2] states 
that a pure entangled bipartite system always shows nonlocality not only for the case 
of dichotomic measurement but also for the case of a measurement with an arbitrary 
number of measurement outcomes. 
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